Squeeze-out of nuclear matter in peripheral heavy-ion collisions 
and momentum-dependent effective interactions. * 
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Abstract 

We perform a systematic study of in-plane and out-of-plane proton and 
neutron flow from nucleus-nucleus collisions within the BUU transport ap- 
proach employing different parameter sets for the mean-field potentials that 



momentum-dependent potential. We find that a simultaneous description of 
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, the experimental data from the BEVALAC and the SIS on both the nucleon 

o 

squeeze-out V2 and the in-plane flow F at beam energies = 0.15 -r- 2 

AGeV requires a mean field with strong momentum dependence, i.e. an effec- 
ts . 

tive Landau mass m ~ 0.68mo at normal nuclear matter density po = 0.17 



fm 3 , where mo = 0.938 GeV is the bare nucleon mass. Some experimen- 



tal data on the squeeze-out require an even stiffer momentum dependence 
(m* ~ 0.62 mo). All systems investigated are found to be compatible with 
m*/m = 0.65 ±0.03. 
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I. INTRODUCTION 



Since many years the study of the various collective flows of nuclear matter is one of 
the main subjects in heavy-ion physics [0-0]. The interest in collective nuclear motion 
under extreme conditions, like high density and/or high temperature, originates from the 
fundamental problem to extract the equation of state (EOS) of nuclear matter. It was shown 
in Refs. fiDHTBf that the in-plane directed flow H] is sensitive to both the incompressibility 
K of the EOS and to the momentum-dependent effective interaction (MDI). Moreover, the 
in-medium nucleon-nucleon (NN) cross section (that is related to the imaginary part of the 



in-medium G-matrix) also influences the transverse directed flow |15| since nucleons acquire 
transverse momentum also by elastic and inelastic collisions. It is necessary, therefore, to 
have at least two additional observables, besides the in-plane flow, in order to determine the 
three most commonly used components of the effective NN interaction adopted in the various 
analyses: the incompressibility K, the stiffness of the MDI and, possibly, the in-medium NN 
cross section. 

The rapidity distribution of nucleons, which is only weakly dependent on the mean- 
field parameters, gives reasonable constraints on the in-medium NN cross section or on 



the nuclear stopping power [|I6| , |T7]| . However, it was not so clear until recently, how to 
possibly extract information on the MDI from heavy-ion collisions. Danielewicz has 
shown within the Boltzmann-Uehling-Uhlenbeck (BUU) transport model, that the elliptic 
flow (or squeeze-out ratio) of protons at midrapidity is very sensitive to the MDI part of 
the nucleon mean field. The reaction Bi+Bi at i?/ a t,=400, 700, 1000 AMeV (measured by 



the KAOS Collaboration |19[) has been analysed in |T8| leading to the conclusion that the 
MDI part of the NN interaction should produce an effective Landau mass m* ~ 0.7mo, since 
momentum independent forces (m* = m ) strongly underpredict the squeeze-out ratio at 
high transverse momenta 



The aim of our present work is: (i) to study the mechanism of the squeeze-out enhance- 
ment by the MDI and (ii) to analyse several independent data sets on the directed and elliptic 
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flow of nucleons to get experimental constraints on the stiffness of the MDI. We concentrate 
on the bombarding energy regime up to 2 AGeV since here multiple meson production or 
"string" degrees of freedom play a minor role. As shown by Sahu et al. [^0| the collective flow 
is strongly influenced by the effective degrees of freedom taken into account in the transport 
approach at high energies. For a detailed investigation of this problem - especially at AGS 



energies - we refer the reader to Ref. |20 



The paper is structured as follows: In Sect. II a brief description of the applied BUU 
model is given. Sect. Ill contains our interpretation of the squeeze-out mechanism which is 
illustrated in Sect. IV by a detailed study of the time evolution of colliding nuclear systems 
within BUU calculations. A comparison to the experimental data on the proton and neutron 
elliptic and directed flows is given in Sect. V, while Sect. VI concludes our work with a 
summary of the results. 

II. THE BUU MODEL 

In our calculations we have applied the recent version of the BUU model described in 



more detail in Ref. [f2T|. The equation of motion for the nucleon phase-space distribution 
function /(r, p, t) includes a propagation in the momentum-dependent scalar mean field 
s(r, p, /) via the single-particle energy 



H mf (r, p, t) = + s(r, p, f)f + . (1) 

Both elastic and inelastic scattering processes are additionally described by the r.h.s. of the 
transport equation 

/ d_ dH mf d _ d(H mf + U Co ui) d \ 

\dt dpi dr dr dpi) ' ^ 

+ coupling terms , (2) 

where Ucoui is the Coulomb potential acting on protons; g = 4 is the spin-isospin degeneracy 
of a nucleon; vu = |vi — V2I and px2 = |pi — P2I are the relative velocity and relative momen- 



turn of colliding nucleons, respectively; d<TNN(pi2)/dfl is the energy dependent differential 
NN scattering cross section as parametrized by Cugnon /, := /(r, Pi,t) (i = 1,2,3,4), 
while P3 and P4 are the final momenta of the scattered nucleons. The coupling terms in (0) 
account for the various inelastic channels: NN ttNN, NN NR - S-wave pion and 
resonance production/absorption, NN — > AKN, NN — > TPKN - strangeness production. 

The scalar potential s(r, p, /) is computed as follows (see [p3fl): 
1) In the local rest frame (l.r.f.) of the nuclear matter element, where the space components 
of the baryonic 4-current = j) vanish, i.e. j = 0, the single-particle energy is calculated 
as 



e/.r./.(r,p) = \]ml + p 2 + U(p,p) . (3) 

The explicit form of the momentum-dependent potential U in (|3]) is taken from Welke et al. 

as 

U(p, p) = A— + B (—\ + — f £%- 3 - / (r ' P ?, 2/A2 • (4) 
Po VPo/ Po J (2tt/i) 3 1 + (p - p') 2 /A 2 

A symmetry energy term is not taken into account in the potential ([|). 
2) The mean field s is calculated from the relation in the l.r.f. 



4r.f. -P 2 ~m . (5) 

3) The scalar potential s is then used in Eq. (JJ) to determine the single-particle energy in 
the computational frame. Steps l)-3) are repeated until selfconsistency is reached, i.e. s is 
not changed within the accuracy of numerics on the level of 10~ 4 . 

The interaction range A of the potential (Q) is chosen as in Ref. 0: A = 1.5^ 0) , where 
pp = Pf(po), Pf(p) = (|tt 2 p) 1 ^ 3 ^. The residual four free parameters A, B, r and C of the 
potential U are determined from the conditions: 
(i) The effective mass at the Fermi surface 

(m*)- 1 = mo 1 + (PfY 1 ^ (0) (6) 

P \p=p F 



must be in the range 0.6mo < m* < tuq (c.f. p5|j26[| ). This condition puts limits on the 
stiffness C of the MDI. 

(ii) The energy per nucleon in nuclear matter has a minimum at p = po and 

(iii) a value of —16 MeV at the minimum: 

^ = 0,1 — -16 MeV , (7) 

Op \p= po P \p=p 

where 

, \ 3 In 2 1 o r+1 

e (p) = -e F p + -A— + ——B^- + 
5 2 p t + 1 Po 

C7 /• ffdp 2 /(Pl)/(P 2 ) /q\ 

Po 7 (2tt^) 3 (2nh) 3 1 + (pi - p 2 ) 2 /A 2 
with = p 2 F /2mo and /(p) = Q(pf — v)- Explicit forms of the momentum integrals in 
(|),(D are given in Ref. @]. 

(iv) The nuclear matter incompressibility 



dp „ ,„„ 



K = 9 Po^ (9) 



is fitted in the range K = 200 380 MeV. 

We have applied in our calculations the following mean fields denoted by H (hard, with- 
out momentum dependence), HM and SM (hard and soft, with a "standard" momentum 
dependence) as well as the new mean-field parametrization which is dubbed HM1 (see be- 
low). The parameters of all interactions are presented in Table 1. The stiffness C of the 



MDI in the interactions HM and SM is close to that of Ref. ||24|| . The set of parameters 
HM1 includes an MDI stiffness that is about 30% larger than in HM and SM, which leads to 
a smaller effective mass at the Fermi momentum for HM1. In Fig. 1 we show the resulting 
energy per nucleon as a function of density and the potential U(p,p) versus momentum for 
several densities. The interactions H, HM and HM1 give practically the same EOS (E/A(p)- 
dependence), but they differ in the MDI part, i.e. in their momentum dependence. This 
gives the possibility to extract an independent information on the MDI stiffness C . On the 
other hand, the interactions HM and SM produce different EOS, but have the same MDI, 
which is relevant for a separate study on the sensitivity to the incompressibility K. 



III. SQUEEZE-OUT MECHANISM 



A squeeze-out of nucleons is observed experimentally in the beam energy range below l-=-2 
AGeV where mean-field effects still play an essential role. It can be interpreted, however, 
to a major extent in terms of a shadowing phenomenon 0, as in the participant-spectator 
picture of a heavy-ion collision the fireball particles - emitted in the reaction plane - are 
rescattered on the spectators. This results intuitively in a squeeze-out ratio Rn '■= (A(90°) + 
A(270°))/(A(0°) + A(180°)) at midrapidity being larger than 1. Here N(<f>) is the azimuthal 
distribution of nucleons with respect to the reaction plane in a given rapidity interval. 
A specific particle emitted from the center-of-mass of a system and moving in transverse 
direction with velocity Vt will be shadowed by the spectator piece, if it reaches the spectator 
during the passage time At = 2R/ (fp' m '7) of the colliding nuclei. Here R is the radius of the 
nuclei (assummed to be equal for projectile and target), Vp" 1 ' is the projectile velocity in the 



center-of-mass system, and 7 = (yl — (v c p m -) 2 ) 1 . A straightforward geometrical estimate 
then leads to the condition v t At/2 > R — b/2, or 

v t > ^r^^r > ( io ) 

where b is the impact parameter. One can see from ( |I0D that at larger impact parameters b 
the lower limit of v t becomes smaller, i.e. more particles will be shadowed and the squeeze- 
out ratio Rn will increase (c.f. later Fig. 12). This is simply related to a larger size (~ b) 
of the spectator pieces in peripheral collisions. It is also evident, that for faster moving 
particles it is easier to fulfill the condition (|II]) and, therefore, the squeeze-out ratio should 
grow with the particle transverse velocity v t (c.f. later Figs. 6,13). 

The squeeze-out phenomenon can also be characterized by a negative elliptic flow v<i (c.f. 
which appears in the Fourier expansion 



N((f)) oc 1 + 2v 1 cos(0) + 2v 2 cos(20) + ... . (11) 
Neglecting thus higher order terms in (|TT|), the following direct relation is obtained: 
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In ( p!l| ) the elliptic flow t> 2 is given as 

* - (4i4) • (13) 

At higher beam energy (~ 10 AGeV), faster moving and Lorentz-contracted spectators 



do not shadow particles any more from the expanding fireball p0 |. This results in Rn < 1 



(or t>2 > 0) as observed in experiment p8 



The simple geometrical picture discussed above will be illustrated by the BUU phase- 
space evolution in the next section. Since the squeeze-out is related to the participant- 
spectator reaction mechanism (c.f. ||29||), we expect a lower beam-energy limit of ~ 0.1 
AGeV for > 1. At even lower beam energies the mean-field effects lead to the formation 
of rotating nuclear systems JJIJ that emit particles predominantly in the reaction plane due 
to centrifugal forces, thus producing R N < 1 (or t> 2 > 0, c.f. later Fig. 7). 

We have to note, however, that the shadowing scenario can not completely explain the 
squeeze-out phenomenon. For instance, a small but statistically significant increase of — Vi 
at the beam energy ~ 0.4 AGeV (see later Fig. 7) with increasing incompressibility K can 
not be understood in the shadowing picture. This is, probably, a manifestation of a focussing 



of high momentum particles by the repulsive mean field as proposed in Ref. [18 



IV. BUU STUDY OF THE SQUEEZE-OUT TIME EVOLUTION 

We study the system Au+Au at the beam energy Ei a b = 0.4 AGeV and impact parameter 
6 = 6 fm. Fig. 2 shows the time evolution of the central baryon density and of the NN 
collision rate for the HM and SM mean fields. In agreement with Fig. 1 (upper left panel), a 
slightly higher density and, as a consequence, a higher collision rate are reached in the case 
of the SM parametrization. However, the various EOS shown in Fig. 1 practically do not 
differ in the density interval < p < 0.27 fm~ 3 probed in the collision. We do not expect, 
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therefore, to get clear constraints on the incompressibility K from our study of peripheral 
collisions and concentrate mainly on the MDI-part of the effective interactions. 

As pointed out in Ref. Jl8[ the pj-dependence of the v-i coefficient is very sensitive to 



the momentum dependence of the nuclear mean field. In order to understand this fact, we 
have performed calculations with the parametrizations H and HM of the mean field. Fig. 3 
shows the average transverse velocity (a) and transverse momentum (b) of the midrapidity 
neutrons moving in- and out- of the reaction plane versus time. We observe that in the 
case of the HM parameter set (vt) reaches a maximum at t ~ 17 fm/c when also the central 
density and the collision rate are close to the maximum (Fig. 2). At this time the system 
shows the most compact configuration in space, i.e. the line connecting the centers of target 
and projectile spectators, is perpendicular to the beam (z-) direction (Fig. 4). However, in 
case of the H mean field the transverse velocity has no maximum (the peaks of (v t ) and (p t ) 
at t ~ 7 fm/c are due to fluctuations, because initially there are no nucleons at midrapidity). 
Moreover, at t < 40 fm/c in the case of the HM mean field the average transverse velocity is 
larger than for the mean field H, while the transverse momentum is practically independent 
on the mean field, even showing a slightly opposite tendency: at t < 40 fm/c the mean 
field H produces a somewhat larger (p t ). This can be qualitatively understood from the 
Hamiltonian equations: 

. _ dH mf _ p | m + s ds 
dp H mf H mf dp ' 
dH mf m + sds 

p = -^- = -^7*- (15) 

where H m f is the single-particle energy (|l|). For the momentum- independent mean field H, 



(nj) reduces to a simple proportionality: v = p/H m f, where H m f ~ 1 GeV. However, for 
the HM interaction the strong momentum-dependent mean field s(r, p, /) is created by the 
density build-up, and the second term in the r.h.s. of ( P|) becomes large (and positive), 
which results in a larger transverse velocity for HM. On the other hand, the force acting on 
a particle depends mostly on the EOS produced by a given mean field or on the pressure 



gradient. Thus, the r.h.s. of (|15|) depends relatively weakly on the momentum dependence 
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of a mean field. 

Later on (v t ) and (p t ) decrease reaching asymptotic values at t ~ 50 fm/c. This decrease 
is due to less energetic products of NN scattering, which populate the midrapidity region at 
later times, while first-chance collision products have larger kinetic energies. Faster moving 
particles (in the case of the HM calculation with respect to H) are shadowed more effectively, 
which causes a larger splitting between in- and out-of-plane transverse velocities for the HM 
calculation at large times. 

In Fig. 5 we show the final azimuthal angle dependence of the midrapidity neutron trans- 
verse velocity (a) and the neutron azimuthal distribution (b). The azimuthal angle modu- 
lation of the transverse velocity is clearly visible, which was first observed experimentally 
||3 1|| . The transverse velocity azimuthal angle dependence and the azimuthal distribution 
of particles have the same shape. The depletion of the particle yield at <fi = 0° and 180° 
is caused by a shadowing of fast moving particles by the spectators; the enhanced yield at 
<p = ±90° correlates with the nonshadowed emission of fast particles. 

Fig. 6 shows the coefficient v 2 vs. transverse velocity (left panels) and vs. transverse 
momentum (right panels) at the two time steps t = 30 fm/c and t = 50 fm/c. Even at 30 
fm/c the particles still feel the influence of the nuclear field and, therefore, the simple relation 
v t — Pt/H m f can be applied at this moment only in the case of the momentum-independent 
mean field H. One can see, that at 30 fm/c the ^(pO-dependence is different for different 
mean fields, while the ^(f^-dependence is very similar for the H and HM calculations. 
This is in agreement with our expectation, that particles moving with the same transverse 
velocity are shadowed in the same way regardless of their transverse momenta. At later 
times, this "transient universality" of vziyt) is destroyed since the influence of the mean 
field will become negligible here and v t = p t j Jm^ + p 2 for both parameter sets. 
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V. COMPARISON WITH EXPERIMENT 



In our comparison to data on squeeze-out observables we now compare in parallel also 
to data on the in-plane flow F, which is defined as the derivative of the transverse in-plane 
momentum component (p x ) with respect to the normalized cm. rapidity y^ = (y/y pr oj)c.m. 
at midrapidity 



F- d ^ ClR) 
F ~^W (16) 

A simultaneous description of both observables, i.e. v% and F, will give a rather strong 
confidence that the dynamics of the participant zone and the interaction of the participating 
nucleons with the spectators are described consistently in our calculation and potentially 
leading to stringent constraints on the EOS and the interactions employed (see also Refs. 
33| , |34|1 where the energy regime below 100 AMeV has been studied in detail). 

Figs. 7-9 show the excitation functions of t>2 and F for protons in case of Au+Au 
collisions in the impact parameter range b = 5 -j- 7 fm corresponding to an estimate given in 
Refs. [p8| , p5| . The data are taken from Refs. |8|j35|-p8| . The time evolution has been followed 
until t max =4:0 fm/c at Ei a b=2 AGeV, while at smaller beam energies t max gradually increases 
reaching 250 fm/c at E lab =25 AMeV. 

The t>2 coefficient (Fig. 7) was calculated for free protons in the cm. rapidity range 
\y\ < 0.1. Free protons were selected by the requirement, that the distance to the closest 
particle within a given parallel ensemble is more than some critical distance d c ~ 3 fm in 
order to separate from protons bound in fragments. As discussed in Sect. Ill, the elliptic flow 
t> 2 vs. Eiab shows a nonmonotonic behaviour: it decreases from positive to negative values 
reaching a minimum at Ei a b ~ 0.4 AGeV, and then it starts to increase moderately again 
with beam energy. In this work, however, we only concentrate on the beam energy domain 
below 2 AGeV and, therefore, the transition from negative to positive v-i at a beam energy 



~ 4 AGeV p8| , p0| is not discussed here. The lower transition energy, i.e. the beam energy 
at which v 2 changes sign, is E TR a — 100 AMeV, in agreement with the FOPI-data p9f. The 
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in-plane flow (Figs. 8,9) increases monotonically in the beam energy region considered. The 
balance energy, i.e. the beam energy at which the flow F changes sign, is Ebal — 50 AMeV 
in our calculation. This value is in between the MSU-data {Ebal — 42 AMeV) and the 
FOPI-data f3j| (E BAL ~ 60 AMeV). 

Now we discuss the effect of the mean field on t> 2 and F. The calculation with a 
momentum-independent mean field H is in a satisfactory agreement with proton in-plane flow 
data below 1 AGeV (Fig. 8), but it fails to reproduce the strong squeeze-out at Ei a b ~ 0.4 
AGeV (Fig. 7). The HM parametrization increases the in-plane flow and squeeze-out (see 
also Fig. 5) and, thus, the in-plane proton flow F is now overpredicted. The lower incom- 
pressibility K in the SM calculation results in a reduced squeeze-out and a reduced flow 
with respect to the HM calculation. Generally, the in-plane proton flow and v-i data are best 
reproduced with the SM mean field. 

In Fig. 9 we compare to FOPI, EOS and recent MSU-47T data on the in-plane flow 
excitation function. The selected data sets are now for a mixture of protons and complex 
fragments and, therefore, they exhibit a stronger flow signal (c.f. Fig. 8, where the data 
are for protons only). A separate study of the fragment flow is out of scope of this work. 
However, in order to understand qualitatively the possible effects from cluster formation, 
we present in Fig. 9 the results of the BUU-SM calculation both for all (solid line) and free 
(dashed line) protons. We find a steeper increase of the in-plane flow for all protons with 
the beam energy from 0.1 to 1 AGeV. The calculation for all protons is in a good agreement 
with data at E^ = 0.4 4- 1 AGeV, but it underpredicts the flow at lower beam energies, 
where the free proton flow F increases earlier with Ei^ and is closer to the data. This can 
be explained by the cluster formation scenario from initially free emitted protons by the 
momentum (rather than coordinate) space coalescence. The balance energy given by the 
SM calculation is E^i ~ 50 (60) AMeV for free (all) protons, which is somewhat larger than 
the value E^i ~ 42 AMeV reported in |38| . But the shape of the in-plane flow excitation 
function at E^ = 25 -i- 60 AMeV is in a remarcable agreement with the MSU-47T data. We 
would like to stress, that the negative in-plane flow in the Au+Au system below Ebal is 
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obtained in our calculations taking into account initial repulsive Rutherford-trajectories (c.f. 
Ref. [ftH , where an opposite result has been obtained). 

In Fig. 10 we present the midrapidity neutron azimuthal distributions with respect to 
the in-plane flow axis for the reactions Nb+Nb and La+La at 400 AMeV in comparison to 
the data from Ref. ||IJ (see Fig. 11 of |K|). In this comparison the BUU events were treated 
as the experimental data |40], i.e. the azimuthal angle of the reaction plane was defined by 



the Q vector: 

Q = E""iS|. ( 17 ) 

V V't I 

where the sum is taken over all charged fragments (i.e. all protons in BUU); vj 7 is the 
transverse velocity of the particle v\ w u = 1 for y(°> > 0.2, w u = for —0.2 < y(°> < 
0.2, w u = —I for < —0.2. Furthermore, we have determined the flow angle by 
diagonalizing the flow tensor Fij in the c.m.s. (see |[41|| ): 

1 v v v u 

Fii = Y. n , ■ (18) 

Finally, a rotation of the cm. coordinate system was performed, first, around the beam 
z-axis to align the x-axis with the vector Q and, second, around the y-axis by the flow angle 
Of (see Refs. [^2[^3j). As shown in Refs. [^2j^3| , the squeeze-out is more pronounced in the 



rotated coordinate system (x',y',z'). The azimuthal distributions in Fig. 10 are presented 
for neutrons with longitudinal momenta —0.1 < (p' z /p'p ro j)c.m. < 0.1. We see from Fig. 10, 
that the data can be reasonably well reproduced adopting the HM1 mean field (see Table 
1), suggesting that the MDI stiffness C in the parameter set HM is not high enough. Our 
conclusion is in line with the results of Ref. [40]], where it was reported that BUU calculations 



with the mean field of (K = 215 MeV, i.e. the SM parametrization) underpredict the 
squeeze-out. 

Fig. 11 shows the in-plane transverse momentum (p x ) of neutrons versus the normalized 
rapidity y^ for the reactions Nb+Nb and La+La at 400 AMeV. We observe that an increased 
MDI stiffness improves the agreement with the neutron in-plane flow data, too. 
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Finally, we have studied the neutron squeeze-out in Au+Au collisions measured by the 



FOPI-LAND Collaboration [ 4"4"| . Fig. 12 shows the impact parameter dependence of the 
squeeze-out ratio Rn for Au+Au collisions at 400 AMeV, which reasonably describes the 
data for the HM calculation (dotted histogram). The HM1 parameter set (solid histogram) 
underpredicts somewhat the squeeze-out ratio in central collisions and overpredicts it in 
semiperipheral (b = 7-^8 fm) collisions. The SM parameter set (dash-dotted histogram) 
and momentum-independent field H (dashed histogram) underpredict the data at large im- 
pact parameters. The calculated i?Ar(6)-dependence has a characteristic bell-shape with a 
maximum at b ~ 7 fm. This agrees with the results of Ref. |[39|| , where the squeeze-out was 



studied for charged particles. The drop of Rn at large impact parameters can be explained 
by a smaller compression of the participant zone and, therefore, smaller transverse velocities 
of particles emitted in peripheral collisions. 

In Fig. 13 the pt-dependence of the ratio Rn is shown for collisions of Au+Au at 400, 



600 and 800 AMeV. In agreement with the previous calculations of Ref. []18| we observe that 
a lower effective mass, i.e. a steeper momentum dependence, increases the squeeze-out ratio 
Rn at large p t (c.f. dotted and solid lines in Fig. 13 d-f). On the other hand, reducing the 
NN cross section decreases the ratio Rn at large p t (c.f. solid and long-dashed lines in Fig. 
13 e). The parameter set H produces a too flat i?Ar(p t )-dependence (c.f. short-dashed line 
in Fig. 13 a). The SM calculation gives the best overall agreement with the data (dash- 
dotted line in Fig. 13 a-c). We see that our BUU calculations reproduce the experimentally 
observed trend of Rn to decrease (at fixed p t ) with collision energy, which was interpreted 



in Ref. [44] as a scaling behaviour of the squeeze-out. 



VI. SUMMARY AND CONCLUSIONS 

The squeeze-out of nuclear matter in heavy-ion collisions at beam energies Ei a b = 0.05-^2 
AGeV has been studied within the BUU approach of Ref. |2T|. As demonstrated, being 
essentially caused by the shadowing of the expanding fireball by cold spectators, the squeeze- 
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out is quite sensitive to the specific velocity of the particles. A mean field with a larger MDI 
stiffness gives a smaller effective mass, and thus a harder velocity spectrum of nucleons. 
In a heavy- ion collision thus faster moving nucleons are shadowed more effectively and the 
squeeze-out ratio Rn is larger for mean fields with stronger momentum dependence. 

The systematic comparison with the experimental data on the excitation function of 
the elliptic flow t>2 |§ favours mean fields with m* ~ 0.68mo, i.e. HM or SM sets. The 
influence of the incompressibility K on the squeeze-out turns out to be quite weak. The 
proton in-plane flow data of the EOS Collaboration at Eiab = 0.25 -j- 1.15 AGeV and the 
fragment flow data of the MSU-47T Collaboration at E\ a b = 25 -=-60 AMeV are consistent with 



the SM mean field. We have found, that a reproduction of the BERKELEY data [40 



on 



the neutron azimuthal distributions and the in-plane flow in Nb+Nb and La+La collisions 
at 400 AMeV requires an enhanced MDI stiffness for p < 1 GeV/c (HM1 parameter set: 
K = 379 MeV and m* = 0.62m ). However, the FOPI-LAND data @ on the neutron 
squeeze-out in Au+Au collisions at 400-^800 AMeV favour the standard MDI stiffness (HM 
or SM parameter sets). 

To conclude, a simultaneous description of the nucleon squeeze-out v 2 and the in-plane 
flow F at the beam energies E\ a b = 0.2 -h 2 AGeV requires a mean field with a strong 
momentum dependence. This corresponds to an effective Landau mass at the Fermi surface 
m*/m = 0.65 ± 0.03 at normal nuclear density. The incompressibility K, however, is less 
well determined since in semiperipheral reactions the average density probed up to 1 AGeV 
is too low. 
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Table 1. Parameter sets for the different mean fields employed in the BUU calculation. 



Not. 


K (MeV) 


m*/mo 


A (MeV) 


B (MeV) 


C (MeV) 


T 


A (fm- 1 ) 


H 


380 


1.00 


-124.3 


71.0 


0.0 


2.00 




HM 


379 


0.68 


-10.0 


38.0 


-63.6 


2.40 


2.13 


HM1 


379 


0.62 


22.2 


29.5 


-82.7 


2.62 


2.13 


SM 


220 


0.68 


-108.6 


136.8 


-63.6 


1.26 


2.13 
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FIGURE CAPTIONS 



Fig. 1 Density dependence of the energy per nucleon in nuclear matter (upper left panel). 
The other three panels show the mean-field potential vs. momentum at the nuclear 
densities 0.5po, Po and 2p . The dashed, dotted, solid and dash-dotted lines correspond 
to the interactions H, HM, HM1 and SM, respectively. 

Fig. 2 Central baryon density (smooth lines, left scale) and NN collision rate (histograms, 
right scale) vs. time for Au+Au at 0.4 AGeV and 6 = 6 fm. The solid (dashed) line 
and histogram correspond the HM (SM) calculation. The horizontal dotted line shows 
the normal nuclear density po for comparison. 

Fig. 3 Average transverse velocity (a) and transverse momentum (b) of neutrons in the cm. 
rapidity range \y\ < 0.1 emitted in the reaction plane (0 = —30° + 30°, 150° + 180°, 
-180° + -150°) and out of the reaction plane (0 = 60° + 120°, - 60° + -120°) for 
Au+Au at 400 AMeV and b = 6 fm. 

Fig. 4 Time evolution of the baryon density in the reaction plane (x, z) for Au+Au at 0.4 
AGeV and b = 6 fm calculated with the HM mean field. Isolines from outer to inner 
correspond to densities in the order 0.05, 0.1, 0.15, 0.20 and 0.25 fur 3 . 

Fig. 5 (a) Azimuthal dependence of the average neutron velocity and (b) neutron azimuthal 
distribution for Au+Au at 0.4 AGeV and 5 = 5 -h 7 fm. Neutrons are selected in the 
cm. rapidity interval \y\ < 0.1. The corresponding elliptic flows are: v 2 = —0.046 (H) 
and v 2 = -0.090 (HM). 

Fig. 6 The elliptic flow v 2 vs. transverse velocity (left panels) and vs. transverse momen- 
tum (right panels) of neutrons in the cm. rapidity interval \y\ < 0.1 for the time steps 
t = 30 fm/c (upper panels) and t — 50 fm/c (lower panels). The calculations have 
been performed with the H (open circles, dashed line) and HM (full circles, solid line) 
parameter sets. The system is Au+Au at 0.4 AGeV and 6 = 6 fm. 



20 



Fig. 7 Excitation function of the elliptic flow v 2 as given by the collection of experimental 
data from various groups || and by the BUU calculations. The BUU curves are shown 
for Au+Au collisions at b = 5^-7 fm for free (unbound) protons with a critical distance 
d c = 3 fm in the cm. rapidity range \y\ < 0.1 (see text). The curves correspond to 
various mean-field parameter sets: short-dashed - H, thick solid - SM, dotted - HM. 
For reference the results of the cascade calculation (without mean field) are shown by 
the thin solid line. 

Fig. 8 Excitation function of the transverse in-plane flow F for Au+Au collisions in com- 



parison to the EOS [35| and E895 [37| proton data. The BUU curves are marked as 
in Fig. 7. The BUU results are impact parameter weighted in the interval 6 = 5 -r- 7 
fm and are obtained for free (unbound) protons with a critical distance d c = 3 fm in 
the normalized cm. rapidity range -0.2 < y<® < 0.3. 

Fig. 9 Excitation function of the transverse in-plane flow F for Au+Au collisions calculated 
with the SM mean field for free protons (d c = 3 fm, dashed line) and all protons (d c = 
fm, solid line) in comparison to the experimental data of the MSU-47T Collaboration 



for Z = 2 fragments [38|] and with the data of the EOS and FOPI Collaborations for 
Z — 1 + 2 fragments |36| . The BUU results are weighted in the impact parameter 
interval 6 = 5^-7 fm. 

Fig. 10 Azimuthal distributions of neutrons with respect to the flow axis (</>') for collisions 
of Nb+Nb and La+La at 0.4 AGeV in comparison to the data (histograms) from 
Ref. |4"0|| . Solid and dashed lines correspond to BUU calculations with HM1 and HM 
mean-field parametrizations, respectively. The BUU results are weighted in the impact 
parameter range b = 1 + 5 fm (Nb+Nb) and b = 1 + 6 fm (La+La) according to the 
estimate b max = 0.5(R P + R t ) in [[4(|. In the calculations we selected free neutrons 
(d c = 3 fm) with longitudinal momenta —0.1 < (p' z /p' P roj)c.m. < 0.1. 

Fig. 11 Average neutron transverse in-plane momentum vs. normalized cm. rapidity for 
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Nb+Nb and La+La collisions at 0.4 AGeV in comparison to the data from |40| . The 
BUU results are weighted in the impact parameter range b = 1 4- 5 fm (Nb+Nb) and 
b = 1 4- 6 fm (La+La); free neutrons with a critical distance d c = 3 fm are selected 
in the BUU calculations. The solid and dashed lines correspond to the HM1 and HM 
parametrizations. The experimental threshold laboratory kin. energy of 55 MeV for 
neutrons [4(J is taken into account which leads to the asymmetric (p x ) (y^) -dependence 



with respect to y^ = 0. 
Fig. 12 The neutron squeeze-out ratio Rn vs. impact parameter b for Au+Au collisions 



at 400 AMeV. The data from Ref. are for neutrons in the rapidity interval 0.4 < 
y/Uproj < 0.6, while the histograms represent the BUU calculations with various mean- 
field parameter sets. Upper panel: short-dashed - H, dash-dotted - SM. Lower panel: 
dotted - HM, solid - HM1. In the BUU calculations free neutrons are selected with 
d c = 3 fm and the angular cuts of the LAND detector are taken into account. The 
calculated azimuthal distributions are shown with respect to the reaction plane given 
by the Q vector (see text). 

Fig. 13 Squeeze-out ratio of neutrons Rn in the rapidity interval 0.4 < y/y pro j < 0.6 as a 
function of the neutron transverse momentum for collisions of Au+Au at 400, 600 and 



800 AMeV. The data are from Ref. p4|| ; the histograms show the BUU calculations 
for different mean fields: short-dashed - H, dotted - HM, solid - HM1, dash-dotted - 
SM, long-dashed - HM1 with reduced (by 30%) NN cross section. The BUU results 
are impact parameter weighted in the region b = 5 + 9 fm, which approximately 
corresponds to the centrality bin E2 for the data. In the calculations the reaction 
plane was determined by the Q vector. 
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